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Multicanonical ensemble simulations for the simulation of first-order 
phase transitions suffer from exponential slowing down. Monte Carlo au- 
tocorrelation times diverge exponentially with free energy barriers Ai^, 
which in L'^ boxes grow as L'^~^ . We exemplify the situation in a study of 
the 2D Ising-model at temperature T /Tc — 0.63 for two different lattice 
manifolds, toroidal lattices and surfaces of cubes. For both geometries 
the effect is caused by discontinuous droplet shape transitions between 
various classical crystal shapes obeying geometrical constraints. We use 
classical droplet theory and numerical simulations to calculate transition 
points and barrier heights. On toroidal lattices we determine finite size 
corrections to the droplet free energy, which are given by a linear combi- 
nation of Gibbs-Thomson corrections, capillary wave fluctuation correc- 
tions, constant terms and logarithmic terms in the droplet volume. Tol- 
man corrections are absent. In addition, we study the finite size effects 
on the condensation phase transition, which occurs in infinite systems 
at the Onsager value of the magnetization. We find that this transition 
is of discontinuous order also. A combination of classical droplet theory 
and Gibbs-Thomson corrections yields a fair description for the transition 
point and for the droplet size discontinuity for large droplets. We also es- 
timate the nucleation barrier that has to be surmounted in the formation 
of the stable droplet at coexistence. 
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1. INTRODUCTION 

First-order phase transitions play an important role in many branches of 
physics ranging from the well known liquid vapor transition to nuclear 
physics protein folding ^ or even to the symmetry breaking in the 
early universe Even for the intensely studied liquid vapor transition 
there are still considerable uncertainties in the calculations of the decay 
rate of metastable states partially due to unknown finite size and finite 
curvature corrections to the free energy and the surface tension of droplets. 

The present study uses the two dimensional Ising-model to investigate 
the influence of such finite size effects, because it is conceptually simple 
and there is a large body of rigorous results, which can be used in the 
comparison to simulation data. The partition function of the Ising-model 
is given by 

conf. 

with the Hamiltonian 

H = Hj - hM := - J2 ^i^j -hJ2^i ' (1-2) 

<i,j> i 

where Hj contains the usual nearest neighbor interaction and the mag- 
netization M couples linearly to a external magnetic field h. We use the 
multicanonical sampling method to study the model in the whole magneti- 
zation interval [—L"^, L"^] at the inverse temperature f3 = 0.7(corresponding 
to T/Tc = 0.63, with Pc = ln(l + V2)/2), which is sufficiently low to pro- 
nounce effects due to first-order phase transitions, but still high enough 
to use a isotropic surface free energy as a good first approximation. Most 
of our efforts will be focused on the magnetic probability distribution 

conf. i 

which up to a normalization factor equals the restricted partition function 
Z{m,L). The distribution Pl{M) was already studied in with 
the aim to understand the dynamics of the decay of a metastable state. 
Our goal in this paper is to produce a quantitative description for high 
resolution Monte-Carlo data of Pl{M) in terms of classical droplet theory, 
including the leading finite-size effects. 

The paper is organized as follows. In section 2 we give a brief review 
of the multicanonical sampling method and present numerical evidence of 
residual exponential slowing down, due to the singularities associated with 
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first order phase transitions. We argue that similar hmitations are present 
in other broad histogram samphng methods. In section 3 we present the 
classical theory for the boundary induced geometric phase transitions in 
toroidal and cube-surface geometries. Section 4 discusses the influence of 
a variety of possible finite size effects. Section 5 gives a detailed account of 
our simulation data for both geometric phase transitions and the droplet 
condensation transition, including the analysis of finite size effects. In 
section 6 we conclude our findings. 



2. MULTICANONICAL METHOD 

Multicanonical (Muca) sampling was invented GSE) iq eliminate the 
exponential slowing down of canonical (Metropolis or heat bath) simula- 
tions near temperature- or field-driven first-order phase transitions. At 
inverse temperatures P > Pc the magnetic probability distribution Pl{M) 
as a function of the magnetization M in the Ising-model has two maxima, 
which we denote by ±M]^"^, separated by a valley, where the probability 
is suppressed by a factor e~^^, due to the additional free energy AF ~ aL 
of the interface present in the two phase region. Already for moderate sys- 
tem sizes canonical simulations are not able to sample these exponentially 
suppressed states and the simulation gets trapped in one of the maxima of 
Pl{M). Similar problems arise in simulations of spin glasses. The Muca 
method remedies this issue by biasing the sampling with a weight factor 
Pl{M)~^, thereby producing fiat histograms G1ISI_ Alternatively one can 
use the inverse of the density of states n{E)~^ for a given internal energy 
E to enhance the sampling of the suppressed states in energy driven first 
order phase transitions. Thus the simulations are performed with an ef- 
fective Hamiltonian H^jf = Hj + f3~^ In Pl{M). Unbiased averages can 
be calculated via 

<A>=^^ — mr^r-^ • (2.1) 

A major practical problem in the application of the method is, that one 
needs a fairly good approximation of Pl{M) to run a efficient Muca sim- 
ulation. The determination of an estimate for Pl{M) by conventional 
importance sampling may consume already large parts of the overall com- 
putation time, especially for systems with a rough energy landscape, like 
spin-glasses. Recursive schemes have been proposed ^ for the cases where 
finite size scaling cannot be used to extrapolate Pl{M) to large system 
size L. Recently several new ideas which tackle this practical problem by 
forcing the simulation into the unfavorable states have been put forward 
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Figure 1: (a) Time series for magnetizations M in Muca simulations at 
/3 = 0.7 on a toroidal 44^ lattice as compared to (b) time series for M on 
a cube-surface SH(20) lattice, which has a similar volume. The statistics 
is 5 ■ 10'' sweeps on the torus and 10^ sweeps on the SH lattice. The 
horizontal lines for the torus indicate the position of shape transitions. 

We now want to focus on the performance of the algorithm, suppos- 
ing we have the exact probability distribution Pl{M) at hand. For this 
purpose we study the system size dependence of the autocorrelation time 
r, defined as the mean time needed by the algorithm to go from the left 
side of the magnetization range at —M^"-^ to the right one at M]^"^ and 
back again. After transforming away the energy barrier due to the inter- 
face we may at best expect a random walk motion of the algorithm in 
the magnetization. A walk between —M and M then takes O(M^) spin 
flips. With M L"^ and the usual definition of Monte-Carlo time in units 
of LF' attempted spin flips (sweeps) we arrive at r ~ for the optimal 
behavior of t(L). The sobering news of Fig. ^ displaying time series 
for the magnetization M in Muca simulations of the 2D Ising-model at 
inverse temperature /3 = 0.7 on a torus and a cube-surface (see Fig. E)) 
- is, that the assumption of random walk like motion in the whole mag- 
netization interval is plainly wrong. For the torus geometry (pbc.) one 
clearly sees that the interval for the magnetization at least is divided into 
three sectors with random walk like behavior, separated by two barriers. 
For the cube-surface, which we abbreviate in the sequel as SH (for surface 
of a hypercube in case ci > 2), the presence of barriers is less obvious. An 
estimate of the magnitude of barriers can be obtained, if r is fitted with 
the form 

T = ArVe^''^^"'^^ (2.2) 

where the maximal surface dVtmaxi which determines the exponential slow- 
ing down, is 2L on a toroidal box and 4(L — 1) on a SH{L) lattice. 
The symbol a denotes the surface tension and V is the volume. We obtain 
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Figure 2: Multicanonical autocorrelation times r for the torus (circles) 
and SH lattices (triangles) in the 2D Ising- model at /5 = 0.7 in units of 
heatbath sweeps as a function of the box volume V. The plot is double 
logarithmic. Curves are explained in the text. 



the values 



from the fits as displayed in Fig. |21 Into the fit enter data, which in 
the figure are displayed with solid symbols. The measured i?-values are 
significantly smaller than values R ^ 1 as expected for non multicanoni- 
cal simulations, but clearly indicate the presence of residual exponential 
slowing down. We show, that the barrier on the torus, which leads to 
a exponential slowing down, is caused by a geometrically induced first- 
order transition from a droplet to a strip domain, whose barrier value 
can be calculated using classical droplet theory. The droplet theory result 
R = 0.1346... agrees within error bars with the measured value of (j2.3p . A 
superficial inspection of the time series for the cube-surface, see Fig. [TJd, 
might lead to the guess, that in this case no barriers are present. We will 
show, that in this case there is a series of three different discontinuous 
transitions in-between phase space regions, where droplets occupy one, 
two, three and four corners on the cube-surface. The barrier values are 



R 



R 



0.121(14) Torus 
0.031(05) SH , 



(2.3) 
(2.4) 
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i?i/2 = 0.02987, i?2/3 = 0.02977 and i?3/4 = 0.03441, which all are smaller 
than corresponding barriers on the torus and, the maximum value i?3/4 
again agrees within error bars with the measured value cited in ()2.4|1 . 



Now, as already pointed out by Leung and Zia one could avoid 
this type of geometrical phase transitions by simulating the Ising-model 
on the surface of a sphere, and actually the cube-surface may be viewed as 
a crude approximation to the sphere. Alas there is no way known to put a 
regular lattice of arbitrary volume on the surface of a sphere, but off lat- 
tice simulation of liquid-gas systems could do the job. Are we then able to 
achieve random walk like Muca dynamics in the whole interval of magneti- 
zations ? Our answer to this question still is no, since there is an additional 
first-order transition present, namely the droplet condensation phase tran- 
sition from a uniform one phase region to the phase separated two phase 
region, where we have yet another essential singularity in the free energy. 
Here again we find a, albeit much smaller barrier for the Muca simulations, 
which is now directly related to the physical nucleation barrier associated 
with the formation of a critical nucleus GMB, Now, what at first sight 
just looks like an embarrassing limitation of the algorithm, at a second 
thought gives us interesting information about the droplet condensation 
phase transition. In the quantitative analysis of this data Gibbs-Thomson 
corrections and other finite size effects will play a prominent role, since 
the volume of the coexisting droplet scales for large system sizes asymp- 
totically as L^/^ and thereby introduces an additional scale. Our findings 
not only are relevant for the algorithmic performance of Muca simulations, 
but most likely limit the performance of other broad histogram methods - 
like the recently introduced Wang and Landau density of states sampling 
method^ - as well, provided the Monte-Carlo covers a portion of phase 
space, which contains one of the mentioned singularities. For purposes 
of illustration we display in Fig. |21 the magnetization density time series 
obtained from a Wang-Landau broad histogram sampling in the magne- 
tization of the 2D Ising-model on a 400^ lattice at f3 = 0.7. The value 
of the Wang-Landau parameter / is / = 1.00007 and the magnetization 
density of states is updated with a heat bath algorithm. The horizontal 
lines in the figure denote the Onsager value uiq and the position rricond 
of the droplet condensation phase transition. The existence of a barrier 
is indicated by the presence of flip flops between two different regions of 
the phase space, below and above rricond- This exponential slowing down 
worsens, if either lattice sizes are increased, or / is tuned to the value 
unity. 
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Figure 3: Time series for magnetization densities m in Wang-Landau 
density of states simulations at /? = 0.7 on a toroidal 400^ lattice and at 
/ = 1.00007. The simulation slows down at the position of the droplet 
condensation phase transition nicond, which we study in the present paper. 



3. CLASSICAL DROPLETS 

Below the critical point the bulk density of the spontaneous magneti- 
zation mo (T) in a infinite system is given by Onsagers solution " , which 
at /? = 0.7 predicts the value mo = 0.99016... for the magnetization den- 
sity. If we restrict the mean density of the magnetization to some value 
in the interval [—mo, mo] the system separates into two phases with mag- 
netizations ±mo(T) divided by an interface. The orientation dependent 
free energy of a interface can be obtained from the spin-spin correlation 
function via a dual transformation ^ . The equilibrium shape of a droplet 
with volume fl of the minority phase, embedded in the majority phase, 
can be obtained by the celebrated Wulff construction ^ . The total free 
energy of the droplet is given by ^ 



J:d = 2Vwn (3.1) 

where 



4 fi^'^o 



= — / dx cosh 
P Jo 



cosh' (2/? J) 
sinh(2/3J) 



— cosh(x) 



(3.2) 
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(a) (b) 

Figure 4: (a) Crystal shapes at the droplet strip transition. With the 
droplet volume fixed at its value L^/vr at coexistence the circular stable 
droplet is deformed into the unstable saddle point lens shaped droplet. 

(b) The percentile increase of the excess length [dQc/s — 2L)/2L as a 
function of of the half base length r/ro in units of tq. At r/ro = 7r/2 one 
reaches the saddle point configuration. The measured data point of ()2.4p 
is plotted with error bars and agrees with the classical droplet result. 



is the volume bounded by the Wulff plot of the orientation dependent 
surface free energy, and <Jo(T) = 2 J + -^Intanh/SJ is the free energy of 
the (1,0) surface, which at /3 = 0.7 has the value ao = 0.89643... . In the 
presence of boundaries different equilibrium shapes can occur in 
the case of a torus with periodic boundary conditions there is a first-order 
phase transition from the droplet shape - as given by the Wulff construc- 
tion - to a strip with surface free energy S5 = 2aQL. The strip is wrapped 
around the torus and the transition point is determined by the condition 
S/) = S5. This was shown rigorously by Shlosman Employing nu- 
merical integration for the integral in Eq. ()3.2|) one can precisely calculate 
the transition point mD/s(T), which at /5 = 0.7 has the value 

mo/s = 0.36974... . (3.3) 

In our Monte-Carlo data the transition is rounded and shifted due to sev- 
eral finite size effects, which will be discussed in detail in the next section. 
For the toroidal and cube-surface geometry we calculate in the rest of 
this section classical equilibrium shapes and saddle point configurations. 
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which determine transition points and barrier heights in-between the dif- 
ferent phases. For sake of simphcity we use a isotropic surface free energy, 
which is still a reasonable approximation at inverse temperature (3 = 0.7 
{T/Tc = 0.63), where the Monte-Carlo simulations are performed. 

Toroidal boxes of linear extent L with periodic boundaries have a vol- 
ume V = and without boundaries the equilibrium droplet shape is just 
a circle of radius r = ^JVl/n, since it minimizes the surface dVl = ^/AnVl for 
a given volume Vt. Leung and Zia ^ argued, that the saddle point config- 
uration at the transition is a lens shaped droplet - see Fig. - formed by 
two arcs with a base length L. At the transition point rriD/s/i^f) = 1 — 2/7r 
this droplet interpolates in-between a spherical droplet of radius tq = L/ti 
and the strip. The arc shape again follows from the minimization of the 
surface at given droplet volume L^/vr. FigureEh displays the excess length 
of the droplet, which is deformed from its spherical shape at r/rg = 1 to 
the lens shaped saddle point configuration with half base length value 
r/ro = 7r/2. The excess surface of the saddle point configuration is ^ 



and we can directly check this value in our Muca simulations, since this 
barrier gives the leading exponential contribution to the autocorrelation 
time and should equal R of ()2.3|) . The value of the classical excess length 
is in fair agreement with the simulation value. Some deviations may come 
from the influence of the anisotropy of the surface free energy, which at 
/3 = 0.7 shifts the location of the transition by a few percents. Cube- 
surface lattices of parameter value L are 2D lattice manifolds denoted 
SH{L), which are identified with the surface of a 3D cube with linear 
extend L — 1. An example is displayed in Fig. where a SH{5) lattice 
with linear extent 4 is represented schematically. A SH{L) lattice has a 
volume V = 6 (L — 2)^ + 12 (L — 2) + 8 and there are eight sites on the 
corners of the cube, which possess three nearest neighbors instead of four. 
A classical droplet configuration at magnetization M = covers four of 
the eight corners and has a minimal surface of length dQ = 4(L — 1) 
separating spin up and spin down domains. A small droplet can lower 
its surface by occupying a corner of the cube. With increasing volume 
it becomes favorable to cover two, three and four corners and similar as 
on the torus there exist shape transitions in-between states, which occupy 
one, two, three and four corners. Due to our choice of a isotropic surface 
free energy all equilibrium droplet shapes consist out of circle segments, 
which are closed around the corners. To systematize our considerations we 
denote the uniquely defined circle segment of volume Q and base length b 



dQp/s - 2L 
2L 



0.1346... 



(3.4) 
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V 



Figure 5: The SH{5) lattice manifold is identified with the surface of a 
5^ cubic box with linear extent 4. 



by S{b, fl) and its arc length by dfl{b, fl). A droplet located at a corner 
has the volume fli = fl[S{b,fl)] — 6^/4 and the surface dfli = dfl{b,fl). 
Note that the base length itself does not contribute to the surface dfl{b, fl) 
of the droplet. The situation is depicted schematically in Fig. Et^, where 
the droplet covers the corner A and the base length has the value b = \/2p 
with p as denoted in the figure. To get the equilibrium shape at constant 
volume fli we minimize dfli with respect to b. This results in 



as compared to dil = V^vrfi for a face centered droplet. Similar consider- 
ations apply to two corner droplets and Fig. IHb displays a droplet, which 
covers two corners A and B. The base length b has the value b = L + 2p, 
again with p as denoted in the figure. For the non-extremal droplet we 
have Q2 = ^[S{b, Q)] — and 8^2 = dQ{b, Q), which again after min- 
imizing with respect to b leads to the surface length 




(3.5) 




(3.6) 



for the two corner equilibrium droplet of volume f22- Equating the surfaces 
fl3.5|) and ()3.6p we find Q1/2 = for the droplet volume where 
the transition from the one corner to the two corner droplet occurs. To 
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(a) (b) 

Figure 6: Generic droplet shapes on SH lattice for (a) one corner and (b) 
two corner droplets. Capital letters denote sites on the corners of the SH 
lattice, small letters denote sites on the edges of the cube and greek letters 
denote sites neither located on edges or corners. Solid lines of the shapes 
contribute to the surface dfl, broken lines do not. Sites labeled with the 
same letter are to be identified. 



calculate the barrier height at the transition we argue, that the saddle 
point configuration between the one corner and the two corner droplet is 
reached, when the site b of Fig. occupies an additional corner: the site 
B of Fig. inb- The volume of this droplet is Q1/2 and its surface has the 
value dni/2 = dQ{b = V2L,3/2L'^), leading to a excess length 

dfli/o — dill 

— = 0.02987 3.7 

4L 

normalized with the length 4L of the "strip" or four corner droplet. Be- 
sides this transition we find two more, from the two corner to the three 
corner droplet and from the three to the four corner droplet. To pa- 
rameterize the possible shapes of a three corner droplet we introduce the 
parameter s = Ad/ L as indicated in Fig. [7| We read of the volume of the 
three corner droplet 

1^3 = ^[S{h = LV(3 + s)2 + (l-s)2, f])] + - 1(1 - s)(3 + s)) (3.8) 

and the surface of the droplet is again the arc length 9^3(6(5), fi). Nu- 
merical minimization of dVL^ with respect to s at fixed ^3 then gives us 
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Figure 7: Generic droplet shape on SH lattice for a three corner droplet. 



the equilibrium shape with three corners inside. Equating d^l2 = d^l^ and 
dQs = dQ^ = 4L we find the volumes ^2/3 and ^3/4, where the transitions 
to a three corner droplet occur. The actual transition values mi/2,m2/3 
and m3/4 are 

^ = ^-, ^ = \ ^ = 0.30237... , (3.9) 

mo 3 mo 3 mo 

where magnetization densities are given in units of the Onsager value mo. 
In Fig.lHlwe display the length of classical droplet surfaces as a function of 
fl/V for one, two, three and four corner droplets and transition points are 
marked by vertical lines. To determine the energy barriers at the 2/3 and 
3/4 transitions we again argue, that for the saddle point configurations 
one needs to deform the 3 corner generic droplet shape of Fig. [7| for fixed 
volume at coexistence in such a way, that either point d of the figure 
collapses onto point A, or that an additional fourth corner is included. 
Both of these shapes are contained in ()3.8j) with s = for the volume f22/3 
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and s = 1 for ^3/4. Using these we find for tlie excess lengtlis 

8^2/3 - dVt2 



4L 



0.02977... (3.10) 



— = 0.03441... 3.11 

4L 

with the largest barrier ()3.11|) being in fair agreement with the value ()2.4|) 
measured in the simulation. 
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Figure 8: Classical droplet surfaces ^fl/V'^^^ on a SH lattice. The curves 
correspond to one, two, three and four corner droplets with vertical lines 
denoting the shape transition points. 



4. FINITE SIZE CORRECTIONS 

In the limit of large box volumes V and large droplets the overall 
magnetization M in the two phase region is created by two areas with 
the Onsager value mo for the magnetization density and opposite sign: 
M = {V — Q)mQ — rriQQ and the droplet volume is a linear function of the 
overall density 

n 1 f m\ , , 

o 1 - — • (4-1) 



V 2 \ mo. 

There exists a variety of finite size effects, which all contribute to the 
droplet free energy and potentially modify relation ()4.H1 for the volume 
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Figure 9: Droplet free energy F^ropiet of ()4.4|) for three m values in vicinity 
and below itlq as a function of the Gibbs-Thomson shift AmcT with < 
ArriGT < ^0 —TTL. With increasing values of m curves move from above to 
below. The saddle points of the droplet free energy are denoted by solid 
triangles, states with vanishing droplet volume correspond to solid circles. 
Curve (b) has a m-value right at the condensation phase transition point. 



of the largest minority droplet on finite boxes. Corrections are due to the 
finite curvature of the droplet surface or due to the restrictions posed on 
fluctuations by the finite size of the droplet or the box, or by the presence 
of several droplets. The literature hosts the following corrections to a 
single droplet ^ 

• Capillary wave corrections of surface string excitations with string 
length L contribute ln(L) and constant terms to the free energy 

with presumably different coefficients in the strip and droplet 

state. 

• Contributions to the free energy due to degeneracies of strip and 
droplet states under translations of the lattice cubic group. On a 
finite 2D system with pbc. a strip has a degeneracy 2L and a droplet 
L^. These entropic terms yield a contribution of order ln(L) to the 
free energy. We expect this type of finite size effect to be present 
only at very low temperatures since the degeneracy is removed by 
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thermal excitations. 

• Tolman corrections to the interface tension due to the small radius 
of curvature of very minute droplets GHMEOl, 

• Gibbs-Thomson corrections ^ describing a shift of the bulk mag- 
netization due to the presence of curved surfaces. 

It is very instructive also, to think of Q of ()4.1|) as an order parameter for 
the formation of an minority droplet of extensive size, which is non-zero in 
the two phase region and vanishes for all states with m > \mo\. There ex- 
ists actually then a phase transition associated with this order parameter, 
at which point the droplet either evaporates into a gas of "small" droplets, 
or the gas condenses. Again for infinite systems the phase transition is 
located exactly at the Onsager mo value, but for finite systems corrections 
are present. Since all the above mentioned finite size effects may conspire 
in the observed shift of the transition it is an interesting and unsettled is- 
sue to study their contributions quantitatively. The droplet condensation 
phase transition has recently been studied with Monte-Carlo simulations 
of Ising-models with fixed magnetization in two ^ and three ^ dimen- 
sions. In these studies the size of the largest minority cluster was measured 
as a function of temperature i.e., along paths perpendicular to those we 
use, measuring at fixed temperature for various values of the magnetiza- 
tion. In d = 2 the authors observed a discontinuous condensation phase 
transition with finite size effects similar to those of Fig. Eland attributed 
them to the Gibbs-Thomson effect. We will quantify this statement with 
the following calculations and with the fits of section 15.31 In three di- 
mensions the finite size analysis for the size of the largest minority cluster 
close to the transition is complicated by the fact, that for small enough 
magnetization there is a infinite percolating minority cluster present al- 
ready in the one phase region ® . The precise influence of this fact at the 
condensation transition deserves further investigation. 

For the description of the condensation phase transition we assume 
that the finite box restricted m partition function of a possibly discontin- 
uous transition can be written as 

Z{m, L) = e-^*°!"='(""'^) + e--^f"'^(""'^) , (4.2) 

up to corrections exponentially small in L with F^uik and F^ropiet being 
suitable free energies in the one phase (bulk) and the two phase (droplet) 
region respectively. In the two phase region we use the classical theory for 
a single extensive droplet, as outlined in the preceding section, amended 
by Gibbs-Thomson (GT) corrections. The bulk phases are described by 
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Ginzburg Landau (GL) theory. The Gibbs-Thomson effect ^ accounts 
for a finite curvature of the surface of the droplet. Microscopically one 
can explain this effect by noting, that the average coordination number of 
a spin at a surface with positive curvature is reduced and thereby the rate 
of detachment into the phase of opposite sign is enhanced and the other 
way round for negative curvature. To lowest order this induces a small 
shift AmcT of the magnetization density of equal absolute value in both 
phases, but of opposite sign. Conservation of overall magnetization then 
leads to a shifted droplet volume 

V 2\ nio mo J ^ ^ 

To calculate AmcT to lowest order we minimize the two phase free energy 

Fdropiet = ctVa^ + csl/Am^r , (4.4) 

where the first part is the surface free energy of a circular droplet of volume 
(size) fl, and the second part is the excess bulk free energy due to the shift 
Aitlgt in an expansion of the Ginzburg-Landau free energy up to second 
order around the bulk value F(±mo), which we choose to be zero for 
convenience. The coefficient C2 = 18.1252318487... at /3 = 0.7 is obtained 
from the low temperature series expansion results of ^ . Minimization of 
()4.4|) results into Attigt oc dQ~^ i.e., a shift of Atjigt proportional to the 
curvature of the droplet, which after reinsertion into the free energy leads 
to a curvature dependent downward correction to the free energy. The 
droplet free energy ()4.4|) is displayed in Fig.|Hlfor three different values of m 
below mo in close vicinity of mo as a function of the Gibbs-Thomson shift 
AmcT, which for each m value in the figure ranges in-between Amor = 
and AmcT = mo—m. With decreasing distance to Onsagers mo the stable 
saddle point solution of situation (a) - as depicted in the figure - turns 
into a metastable one (b), which then turns unstable in situation (c), very 
close to mo- 

In finite systems the otherwise stable classical minority droplet of the 
two phase region becomes metastable at the condensation phase transition 
point mcond{L), where the system performs a finite size rounded transition 
into the one phase region. The condensation point shift due to the finite 
system size 

Amcond{L) =mo- m^ondiL) , (4.5) 

corresponds to situation (b) as depicted in Fig. IHl and can be calculated 
by equating the free energy 



Fbuik = VciAml^^^ 



(4.6) 
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of the bulk state - the sohd circles in Fig. El- to the saddle point free energy 
()4.4|) in the two phase region: the solid triangles in Fig. IHl Note that the 
surface term in ()4.4|1 depends also on /S.mcond- One finds the finite system 
condensation phase transition point 

Am,o„d(L) = y4eo„dL-2/3 (4.7) 

on toroidal boxes. The exponent —2/3 for the finite-size behavior of 
the condensation transition was already found before in the context of 
metastable decay ^ but the inclusion of Gibbs Thomson corrections 
lead to a different value of Acond- The coefficient Acond at /? = 0.7 has the 
value 

A,ond = mor]—^ = 0.23697... , (4.8) 
where parameter values rj and a are given by 

V = [4^]^^^ = 0.20102... (4.9) 
a = y/W/n = 0.90358... . (4.10) 
The droplet size Qcxc{L) at the condensation phase transition is 

f]^^^(L) = = 0.07977... L^/^ ^ (4^1^) 

3 mo 

and one can also calculate the nucleation barrier Bnucu which at the 
condensation phase transition corresponds to the maximal excess droplet 
free energy in-between states at the saddle point and states at Attigt = 
mo — m. One finds 

nud 0.174038... (4.12) 



in units of the coexisting droplets surface free energy and Bnuci / Fdropiet = 
0.154701 in units of the total free energy of the coexisting droplet. The 
Gibbs-Thomson corrected droplet looses its stability at the point 

Am^r^stabieiL) = tjiovQ'^' d)'^' L~^'^ (4.13) 

o 

= 0.18808... L-^'^ , (4.14) 

which inside of the metastable region terminates the droplet metastable 
branch. Condensation phase transition point scaling with the power law 
can also be proven with rigorous methods (see ^ for a recent 
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overview). Furthermore one finds at the condensation phase transition 
point the hnear relation 

AmoTiL) = ^ArricondiL) , (4.15) 

which expresses the Gibbs- Thomson shift in units of the finite size location 
of the phase transition point. This linear relation is especially interesting 
as the ratio 

_ VAmcondiL)/2mo 

^ ~ nUL) ^ ' 

is predicted to have the value Q = 3/2. If one neglects Gibbs-Thomson 
corrections altogether one finds Q — 1. The value of Q does not depend 
on the particular values for the surface free energy, nor on the parameter 
value C2 (which all factor out in the calculation). A stringent test on 
the presence of Gibbs-Thomson corrections thus can be devised, if in the 
Monte-Carlo simulation one measures both, the coexisting droplets size 
^cxc{L) and location of the condensation phase transition ArricondiL) . 



5. SIMULATION RESULTS 

We have simulated the 2D Ising-model at /3 = 0.7 {ksT = 1.428 . . . ) 
on toroidal lattices and on cube-surface SH{L) lattices. In one set of 
Muca simulations we cover all states from the strip to the bulk phase 
and study L = 12 up to L = 44 toroidal boxes and cube-surface 
lattices of parameter values L = 4 up to L = 26. In another set of 
simulations the condensation point phase transition is studied on toroidal 
lattices of sizes 40^ up to 400^, again also with the use of Muca simulations 
covering however a smaller m-interval in vicinity of the condensation phase 
transition. 



5.1. SHAPE TRANSITIONS 

In order to study the droplet shape transitions on the torus and on 
SH lattices in detail we consider the constraint magnetization ensemble of 
the Ising-model. Its partition function Z{m, L) is 

Z(m,L) = ^e-^^<5(m-i^.,) (5.1) 

conf. i 

and we determine the expectation values of observables at fixed magne- 
tization m. Muca ensemble simulations are perfectly suited for the eval- 
uation of restricted expectation values. All one has to do is to average 
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over operator values at m, if the magnetization m is visited in the Muca 
simulation. As we are dealing with droplets it is sensible to introduce the 
m-dependent functions 



dnoim) = J27rl^(l ) (5.2) 

V mo 

fio(m) = ^(1-—) , (5.3) 
2 mo 

which for large system volume V denote the classical surface and volume of 
a circular droplet at fixed overall magnetization m. For each configuration 
at m we evaluate connectivity components i.e., clusters of spins. Nearest 
neighbor spins belong to the same cluster, if they have the same value. The 
size Q - the number of sites - is determined for each cluster and clusters 
are sorted with respect to their size. The second largest cluster defines 
the object of interest and corresponds to the minority phase droplet. Its 
expectation value in the magnetization bin m is denoted <Q> (m) , the 
"volume" of the minority phase droplet. It should be noted that <Q> (m) 
not necessarily has to agree with QQ{m), if finite size effects are present. 
For the toroidal lattice geometry (pbc.) and for the minority phase droplet 
a rectangular bounding box with linear sizes Li and L2 is determined in 
such a way, that the droplet exactly fits into the box. From the geometric 
numbers Li and L2 two geometric order parameters are formed, namely 

Oi(m) = 5(L-max(Li,L2)) (5.4) 

max(Li,L2) , . 

02{m) = . (5.5) 

Values for the operator Oi(m) are zero and unity. At the position mo/s 
of the droplet strip transition we expect to find a finite size rounded jump 
of the expectation value <Oi(m) > from value zero to unity. The order 
parameter 02(m) on the other hand has several discrete values in the in- 
terval [0, 1]. One can construct a probabihty distribution function P{P2) 
for the occurrence of values O2 in the restricted partition function ()5.1|) at 
m. If our classical arguments on the nature of the droplet to strip tran- 
sition are correct, then one expects to find a double peaked distribution 
function P{02) in the vicinity of mn/s with a barrier, which is related to 
the excess length. Susceptibilities Ci(m) and C2(m) are defined by 

Ci(m) = L <(Oi(m)- <0i(m)>)2> (5.6) 

and similar for C2(m) with O2, where < > again denotes the expectation 
value at given m. These susceptibilities show a finite size rounded peak 
at the shape transition defining finite size shifted m£)/5 (L)-values. 
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4<0,>(m) C,(m) <03>(N^,m) 



12 - 




Figure 10: (a): Droplet to strip transition geometric order parameter 
4 < Oi(m) > and susceptibility Ci(m) for the torus with pbc. and (b): 
geometric order parameter <03(m, Nc) > on a SH[2Q) lattice with Nc = 
0, 1, 2, 3 and Nc = 4. The 5 different A^c-dependent curves have maxima, 
which from the right to the left correspond to values A^^c = 0, 1, 2, 3 and to 
Nc = 4. Vertical lines denote classical shape transition points. 

For the minority droplet on SH lattices we define a geometric order 
parameter sensitive to the number Uc of corners occupied by the droplet. 
At magnetization m we count the number Uc and define 

0,{m,Nc) = 6{Nc-nc) , (5.7) 

which at given m receives unity contributions only if Nc corners are oc- 
cupied. We expect e.g., that expectation values <03(m, 2) > yield non- 
vanishing values < 03(m, 2) >^ 1 only, if the magnetization density m 
lies in-between the 1/2 and 2/3 shape transition magnetization density 
values. Data for the geometric order parameters <0i> (m) on toroidal 
lattices and data for < O^lm, Nc) > with iVc = 0, 4 on a SH{26) lattice 
are displayed in Figs. ITUk andllUb. The maximum positions of the sus- 
ceptibility in Fig. [TUk on toroidal lattices define finite volume estimates 
of the magnetization density mu/s at coexistence. Finite size correction 
are not particularly small and thus a careful infinite volume extrapola- 
tion is needed. For every pair of values L and m£)/s{L) we calculate 
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Figure 11: (a): Fit to m^)/^ on the torus. The function OVLq/L is eval- 
uated at the measured finite size shifted values of niD/siL) and plotted 
as a function of L. The line corresponds to the exact result, while the 
solid triangle corresponds to the infinite volume extrapolation. The curve 
corresponds to the fit as explained in the text, (b): Fit to the three shape 
transition points on SH{L) lattices. Here the function y/O^JV is evalu- 
ated at the measured finite size shifted magnetization density values and 
plotted as a function of L. Again horizontal lines correspond to droplet 
calculations, while curves correspond to fits explained in the text. Solid 
symbols correspond to the infinite volume extrapolation. 



dQQ{m£,/s{L))/ L and with the use of the extrapolation 



dno{mn/s{L))/L = /27r(l - ^) + Ao/s^^ (5.8) 

mu/s is determined. We mention, that ln(L)/L finite size corrections 
are predicted by capillary wave fiuctuation corrections to the droplet free 
energy ^nd by entropic terms. The fit in accord with ()5.8p has a 

X^o/ value 0.32 and results into a transition point at mo/s = 0.376(9), 
which within error bars coincides with the exact result ()3.3j) . The data 
and the fit are displayed in Fig. ITTk . A similar analysis based on C2{m) 
data yields the determination mo/s = 0.360(10), which again agrees with 
the exact result. Similar as in torus case we also calculate susceptibilities 
of the O3 order parameter on SH{L) lattices and determine finite size 
shifted shape transition points. The results are displayed in Fig. ITTb as 
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Figure 12: (a) Geometric order parameter O2 distribution on the torus 
in vicinity of ttid/s and (b) barrier values as a function of 1/L. The 
horizontal line denotes the classical droplet theory result, the solid triangle 
the infinite volume extrapolation. 



a function of L on SH{L) lattices. In this case three different i = 1,2,3 
infinite volume transition points are determined via the extrapolation 

V — V — - y - ^ ^'^'^'L ■ ^^-^^ 

We did not include a ln(L)/L term in this Ansatz, since capillary wave 
corrections are the same on both sides of the transition and entropic 
terms should be absent since the droplets center of mass is energetically 
pinned. We obtain the values mi/2 = 0.661(2) mo, ^2/3 = 0.342(3) mo 
and m3/4 = 0.305(3) mo, which within a systematic one percent relative 
error all agree with the values of ()3.9j) . The small discrepancy is presum- 
ably caused by the isotropic surface free energy of our (approximative) 
droplet calculation for SH lattices. The susceptibility C2(m) on toroidal 
lattices has its maximum value at positions mD/s{L) and we can deter- 
mine the probability distribution P{02) of the order parameter O2 there. 
The data are displayed in Fig. IT^ for L = 20, 28, 36 and L = 44 lattices. 
One observes clear double peaks with the peak for the strip concentrated 
in the single point at O2 = 1. In-between one finds states, which with 
increasing lattice size become less an less probable. These suppressed 
states are the saddle point crystal shapes of Fig. HI A split of phase space 
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into two disconnected regions with a free energy barrier in-between proves 
the discontinuous nature of the droplet strip transition and one can de- 
termine the barrier height. The distribution functions have two maxima 
values Pmax,i and Pmax,2 and a minimum value Pmin- We form 

^Fdroplet = -^T^{Pmax,lPmax,2/ Pmin) ' (5.10) 

which we measure in units of AFst^ip = ^Lcxq and the result is displayed in 
Fig. as a function of A linear fit in (the curve in the figure) 
results into the value 

AFdroplet/AFstr^p = 0.133(6) (5.11) 

at a Xdof value 0.40 for the fit. The free energy barrier measured in 
exponential slowing down of Muca simulation (j2.3|) agrees with the free 
energy barrier measured in the suppression of saddle point crystal shapes 
()5.1H) and with the result of classical droplet theory ()3.4|) . 

5.2. MORE DROPLET FREE ENERGY CORRECTIONS 

A direct consequence of a discontinuous behavior of the restricted 
partition function Z{m,L) at the droplet strip shape transition rriD/s on 
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the torus is again the vahdity of a superposition Ansatz ^ 

Z{m, L) = e-^=*"f("^'^) + e--^'^-f'='(™'^) , (5.12) 

which we use to determine free energy corrections. If we define the quan- 
tity 

X = = - 1 (5.13) 

we observe, that classical bulk and surface free energy contributions pro- 
portional to and L cancel in Fstripijn = 0, L) — Fdropietijn = rnu/s, L). 
Note that the classical surface contributions to the free energy of a strip 
state at m = equal those at m = mjjjs and in addition classical surface 
portions of Fgtrip equal those of F^ropiet at m = rriD/s- We obtain the 
representation 

— ^Gstrip{m=0,L)-Gdropletijn=mo/g,L) j^^-j 

where the L-dependence of the function G at most is of the order o(L) 
and possibly contains the whole set of finite size corrections to strip and 
droplet states. 

One of the quantities, which in the Muca simulation is determined, 
is the magnetization probability Pl{M). It counts the probability to find 
magnetization M in the unconstrained Ising model and is proportional to 
the constraint partition function Pl{M = mV) oc Z{m,L). The value 
Pl{M = 0) is easily measured at the magnetization bin M = 0, while 
Z{m = mD/s,L) with mo/s = 0.36974... is interpolated from data at 
magnetization bins closest to the value Vuid/s- Fig. [121 displays the data 
for X as a function of L. Right at the crystal shape transition and with 
increasing lattice size droplet states relative to strips turn out to be more 
and more probable. For large systems the data can be fitted with a power 
law in L (the solid symbols in the figure) 

X = AxL"" (5.15) 

and at a x^^j value of 0.12 we obtain with Ax = 1-4(4) an exponent value 
of a = 0.44(8). Such finite size effects do not have their origin in Gibbs- 
Thomson corrections but either are generated by capillary wave fluctua- 
tion corrections - or by the count of strip and droplet states with respect 
to translations. At zero temperature each translational degree of freedom 
reduces the free energy by a term — ln(L) and for finite temperature, each 
fluctuating surface contributes +| ln(L) via the capillary wave expansion. 
Our theoretical prediction for a therefore is a = — 1 + 2(|) + 2 — ^ = 3/2, 
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Figure 14: Finite size rounding of the discrete derivative (AmlnZ)/L in 
vicinity of the crystal shape transition (a) and in the droplet phase (b). 
For reasons of clarity we only present data for two lattice sizes, 20^ and 
40^. The curves labeled "c" in the figures exactly match the measured 
data (circles) and correspond to theoretical predictions for droplet and 
strip free energies. The data denoted by crosses in (b) are results from 
a simulation with a modified partition function (one spin is fixed to — 1) 
and agree with curves labeled "b". These are based on the droplet free 
energy without logarithmic phase space factor in the droplet volume. 



if translational degeneracies are counted -ora = 2(|) — | = l/2, if tem- 
perature lifts degeneracies. The measured value is consistent with one half 
and thus temperature is too high for a "naive" zero temperature count of 
degeneracies. Note, that two independent surface strings were assumed 
for the strip and one for the droplet. To probe free energies further, we 
introduce the discrete m partition function derivative 

^mlnZ:= ^^-j^^ , (5.16) 

which for small Aj\/ is proportional to ^^^\t,v- Thus AmlnZ/L^ can 
be interpreted as a magnetic field or as a chemical potential ^ in the 
lattice gas interpretation. We have chosen a suitable value = 20 and 
the measured Monte Carlo data are displayed in Figs. IT^ and 114b for m- 
values in vicinity of m = mjj/s {^) ci^nd for m values in the droplet phase 
in (b). Similar data for the 2D Ising model at temperature T = 0.8Tc, 
presented in Fig. 10 of also showed the presence of the strongly rounded 
stripe droplet transition, without attempting a quantitative fit of the data. 
For reasons of clarity we only present results for two lattices, 20^ and 40^ 
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in size. The data, the circles in the figures, can be compared to several 
finite size rounding predictions (the curves in the figures), which are based 
on the superposition Ansatz ()5.12|1 . The infinite volume prediction, which 
is labeled "a", jumps at rriD/s from the value zero (strip phase) to a finite 
value (gap) and in the droplet phase follows the classical result flH.l|l . The 
curves labeled "b" in the figure correspond to finite size rounding in accord 
with ()5.12|) and with free energy functions 



including classical terms as well as capillary wave fiuctuations corrections. 
As can be noted, the presence of further correction terms is suggested. 

The translational invariance of a droplet of — spins fioating in a back- 
ground of + spins can easily be broken, if one spin of the partition function 
Z droplet is fixed to the value —1, which in consequence lowers the partition 
function by a factor / : Zfi^^d spin = fZdropiet with / < 1. At zero temper- 
ature one finds / = fl/V and thus a "microcanonical droplet phase space 
volume" correction to the droplet free energy of the form 



quite similar to the microcanonical phase space volume of gases, is pre- 
dicted. Adding this term to the droplet free energy in the droplet phase 
moves curves labeled "b" to curves labeled "c" in Fig. H^b, which exactly 
reproduce the data. One can also do a simulation in the modified theory, 
which differs from the original one by the fixation of one arbitrary single 
spin to the value —1, which never is updated. The crosses of Fig. 114b 
correspond to data from such a simulation and as can be seen: they come 
to lie on the curves labeled "b" , the free energy form without logarithmic 
phase space factor in the droplet volume. 

It is interesting to ask the non-trivial question, whether similar phase 
space corrections contribute to the strip free energy. Strips are sepa- 
rated from droplets by barriers and thus the full occupancy of phase space 
(with logarithmic phase space factor in the droplet volume) may never be 
reached. A situation like this can be termed: spontaneous breaking of 
translational invariance for strip states, which according to our findings 
actually is realized in the f3 = 0.7 2D Ising model. The curves labeled (c) 
of Fig. have been calculated with the droplet free energy Fdropiet and 
with 



Fstrip = 2o-oL + ln(L) 



(5.17) 
(5.18) 



Fdropiet = aVAjdi + C2V Aml,rp + -ln{dQ) 




(5.19) 



F strip Fgi'^ip In (2) 



(5.20) 
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Figure 15: (a) Logarithmized magnetization density probability distribu- 
tion data as a function of m at /? = 0.5(T/Tc = 0.88) on a 20^ box and 
finite size rounding fit as explained in the text. Data and the fitted curve 
lie on top of each other. The fitted infinite volume interface tension ctq is 
displayed in (b) as a function of The curve in (b) corresponds to the 
exact result. 



and coincide exactly with the data. The strip free energy lacks logarith- 
mic phase space factors, which are inconsistent with the measurement and 
we have included a constant correction term. Each droplet at the crys- 
tal phase transition tunnels into either one of two possible different strip 
configurations and thus the strip free energy is lowered by a term — ln(2). 

Droplet and strip free energies of Eqs. ()5.19p and ()5.20|) can be used 
to determine the interface tension ctq from the finite size rounding of the 
free energy at the crystal shape transition. If for practical purposes we 
approximate the droplet free energy by 

Fdropiet ~ adVLQ{m) + + ln(— ) (5.21) 

we obtain with Eq. ()5.12|) and with lnZ(m, L) = InPL^m) + const at fixed 
mo a four parameter representation of the free energy as a function of 
the parameters ctq, a, ugt and const, which easily can be fitted. A typical 
set of data for the logarithmized free energy at /? = 0.5 on a 20^ box is 
displayed in Fig. IT^ and the fit for values m < 0.8mo (droplet and strip 
phases) lies right on top of the data. Figure ITSb displays the interface 
tension ctq as determined from rather small boxes as a function of {3 in 
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Figure 16: Logarithmized magnetization density probability distribution 
functions Piifn) on toroidal boxes at (3 = 0.7 for the 2D Ising model. 
The maxima are normalized to values unity in Pilm). One observes al- 
most Gaussian fluctuations in vicinity of mo. Once one approaches the two 
phase separated phase space region, a finite size rounded cusp structure 
appears, which corresponds to the condensation phase transition. 



comparison to the exact result, the curve in the figure. The agreement is 
excellent demonstrating, that finite size corrections of the free energy are 
faithfully represented. 

5.3. CONDENSATION PHASE TRANSITION 

The phase transition of an extensive minority droplet in the phase 
separated phase space region into a gas of small droplets in the bulk phase 
is studied on toroidal lattices. In order to provide an overview over the 
data we display in Fig. the (logarithmized) probability distribution 
PL(m) of the magnetization density - the constraint partition function 
Z{m, L) - as a function of m in vicinity of Onsagers magnetization density 
value tuq. A finite size rounded cusp structure is visible for values of m 
slightly below mo and corresponds to the position of the condensation 
phase transition mcond{L). The data are obtained on 40^ up to 400^ boxes 
with the help of Muca simulations, which quite similar to simulations with 
the Wang-Landau algorithm (see section |2)) also suffer from exponential 
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Figure 17: Expectation values <Q> {m)/L'^ for the minority droplet size 
density in vicinity of the condensation point phase transition for 40^ up 
to 400^ boxes. The infinite volume prediction (1 — m/mo)/2 corresponds 
to the dashed diagonal straight line. 



down and in fact: it would be quite time consuming to obtain data of 
comparable statistical quality for a 1000^ box. 

The expectation value for the size of the minority droplet <Q> (m) 
is displayed in Fig. El From the peak positions of the Q susceptibility 
(not displayed in a figure) we determine finite size shifted condensation 
point magnetization values mcond{L) and calculate the shift due to the 
finite system size AnicondiL) = mo — mcond{L). The shift is displayed in 
Fig. El in a double logarithmic scale as a function of L. Numerical values 
are contained in table ^ The condensation point shift of Fig. El contains 
a straight line, which corresponds to the theoretical prediction of classical 
droplet theory with Gibbs- Thomson corrections ()4.7|) . While finite size 
corrections to the theory on small boxes are large, it appears perfectly 
possible, that the data (the solid circles in the figure) approach the theo- 
retical prediction on large boxes i.e., for large droplets. Additional finite 
size corrections are caused by sub-leading free energy contributions and 
quite similar to the discussion of section 15.21 we compare in Fig. El data 
for the discrete derivative A^lnZ on 200^ and 400^ boxes with finite size 
rounding theory, as predicted by the superposition Ansatz (j4.2j) . Simi- 
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Figure 18: Condensation point phase transition transition points at (3 = 
0.7 in the 2D Ising model. The straight hne corresponds to the theoretical 
prediction, including Gibbs-Thomson corrections. The curve corresponds 
to a fit to Eq. ()5.24|) as explained in the text. 



lar Monte-Carlo data for systems of smaller linear extent were shown in 
Fig. 10 of for the Ising model in D=2 and in Figs. 6 and 7 of ® for 
D=3. The onset of the same singularity as displayed in Fig. ^Jis clearly 
visible in those data, despite of the presence of strong finite-size rounding. 
The Ansatz now uses the precisely known droplet free energy of (|5.19p and 
the bulk free energy of (|4.6|) . The dashed lines of the figure do correspond 
to finite size rounding predictions. It is evident, that the second order 
expansion of the Ginzburg-Landau free energy is not precise enough. The 
calculation of coefficients Cn 



of an polynomial expansion for the constraint effective potential Ueff{m) 



oo 




(5.22) 



n=2 




(5.23) 



in higher powers of m — mo, could provide a better bulk free energy form. 
We will not however pursue this issue here. 
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Figure 19: Finite size rounding of data for the discrete derivative 
( A^lnZ) /L as defined in Eq. (j5.16|) for m values in vicinity of the conden- 
sation phase transition on two lattice sizes, 200^ and 400^. The dashed 
curves correspond to finite size rounding predictions. Discrepancies are 
observed for bulk states slightly above the finite system phase transitions. 
The effect is due to the second order Ginzburg-Landau free energy expan- 
sion. 

We perform a three parameter fit to AnicondiL) data and parameterize 
additional finite size corrections through a single power correction with a 
free exponent value /? 

Am,ond{L) = A,ondL-^^^ + BL'" . (5.24) 

The fit corresponds to the curve in Fig. [THland at a Xdof value of 0.25 we 
obtain the fit parameters B = 23(9), /? = 1.93(11) and Acond = 0.237(3). 
The fitted value for Acond perfectly agrees with the theoretical prediction 
Acond = 0.23697... fl4.8|) . which at the condensation phase transition yields 
support for the validity of classical droplet theory and the presence of 
Gibbs-Thomson corrections for large droplets. A somewhat more stringent 
test of the presence of Gibbs-Thomson corrections can be devised, if at 
the condensation phase transition the ratio Q of ()4.16|) is calculated. For 
large droplets one expects, that Q equals Q = 3/2. For magnetization 
values at - or close to - the condensation transition, we determine the 
probability distribution Pl{^) for the occurrence of a minority droplet 
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0.0976(8) 


45.3(4) 


1.555(25) 


0.013(2) 


200 


0.007762(50) 


0.0878(5) 


102.7(6) 


1.526(12) 


0.045(4) 


400 


0.004575(12) 


0.0825(9) 


243.4(6) 


1.518(16) 


0.085(4) 


oo 


0. 


0.07977... 


oo 


1.5 


0.17404... 



Table 1: Measured observables for the 2D Ising model droplet conden- 
sation phase transition at /? = 0.7. The lowest line contains the infi- 
nite volume theoretical predictions of classical droplet theory with Gibbs- 
Thomson corrections. 



of size Vt in the magnetization constraint partition function at m and 
tune the magnetization values to the point, where finite system double 
peaked probability distributions have equal height. The result Pl^S^ is 
displayed in Fig.EOlas a function of ^jL^I^ for boxes 100^, 200^ and 400^. 
Our numerical simulation demonstrates the existence of double peaks in 
Pi^iyi). There also exists a barrier for states in-between, which increases 
with increasing droplet size. These findings provide numerical evidence 
for the fact, that the droplet condensation phase transition of the 2D Ising 
model is of discontinuous nature. Upon fitting Gaussians to the right hand 
side peaks of Fig. 1201 we determine values for the coexisting droplets size 
^cxci which are contained in rows three and four of table Q Values for 
the ratio Q are given in the fifth row of table H On the 400^ box we 
find Q = 1.518(16), which indeed is very close to Q = 3/2 and thus again 
provides solid support for the asymptotic correctness of classical droplet 
theory and the presence of Gibbs-Thomson corrections. 

A very interesting quantity is the nucleation barrier Bnuci, which at 
the condensation phase transition determines the suppression of states in- 
between the equilibrium minority droplet - and the gas of small droplets 
in the bulk phase. At equal height of the probability distribution Bnud 
can be estimated through 

Bnud ~ In , (5.25) 

-'mm 

where values Pmax and Pmin denote maxima and minima of Pl{Q). The 
last row of table Q contains Bnud values, which are given in units of the 
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p,(n) 




Figure 20: Probability distributions Pl(^) of the minority droplet size 
at the condensation phase transition. The vertical line corresponds to the 
"gap" prediction of classical droplet theory with Gibbs-Thomson correc- 
tions. 



droplet surface free energy ay^AirQcxc- The data exhibit large finite size 
effects and most likely one needs much larger droplets in order to draw def- 
inite conclusions. On the 400^ box we find a barrier, which is smaller by a 
factor of 1/2 than the theoretical prediction ()4.12|) . Chances thus are high, 
that the condensation of large droplets proceeds within Gibbs-Thomson 
corrected classical droplet theory through the formation of intermediate 
saddle point configurations, which are classical in nature also. 
The grand-canonical droplet model of GHSJ predicts 



crit 



{d-l)a 
2dmoh 



d> 2 



(5.26) 



for the volume ^crit of the instable critical droplet as a function of the 
ordering field h. Taking into account that the the surface tension a = 
v^ivrcr of EMU differs from ours by a factor v^ivr one obtains 



770" 



crit 



(5.27) 



in d = 2. Using Eq. 1)4.111) and Am = h/2c2 we compare Qcrit with the 



34 



Neuhaus and Hager 




Figure 21: Probability distribution P{E) for the energy in the 2D q = 20 
Potts model on a 70^ box. The conjectured positions of condensation 
phase transitions are indicated by finite size rounded cusps. 



volume of the equilibrium droplet at the condensation phase transition 

which has the same field dependence as fl5.27|) . ^Icxc is a factor nine times 
larger than Qcrit- As one can see in Fig. 1201 this is in accordance with our 
simulations. The droplet volume ratios of the stable droplet at coexistence 
over the critical droplet at the saddle point approximately have values 2, 3 
and 4 for the considered lattice sizes. It is very well possible, that these 
volume ratios approach the predicted value 9 in the infinite volume limit. 

Finally we want to discuss how the free energy = In Pl{M) in 
general and the observed transitions in particular affect the dynamics of 
the decay of a metastable state. As it stands F{L) describes the equi- 
librium state for conserved order parameter m. It tells us for example 
whether adatoms with a given density Am on a surface of size form 
a large single crystal island or a homogenous adatom gas. For systems 
with conserved order parameter, like binary alloys the typical experimen- 
tal setting is a rapid quench from a uniform state at high temperature 
into the two phase region, followed by nucleation of droplets or spinodal 
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decomposition. In magnetic switching experiments ^ a system with a 
non-conserved order parameter is driven from one equihbrium bulk state 
to the other by an external magnetic field h and pending on the field 
strength one observes two different tunneling regimes: In the stochastic 
regime at small fields tunneling proceeds via the nucleation of a single 
droplet and as the process is rare one observes large fluctuations in tun- 
neling times, if the experiment is repeated many times. At large fields 
multi-droplet states dominate the intermediate configurations of the early 
tunneling process. The multi-droplet nature of the tunneling in this case 
reduces fiuctuations of tunneling times. The crossover region in-between 
both regimes was called the dynamic spinodal 

The decay of supersaturated bulk states in any case proceeds via 
the formation of a single droplet - or via an ensemble of several critical 
droplets. Let us consider a subsystem of supersaturated matter of linear 
extent ^ and volume at fixed finite h, which is embedded into a larger 
system of finite linear size L and volume L^. The properties of critical 
droplets in vicinity of the condensation phase transition determine cer- 
tain aspects of the tunneling processes, which can be activated within the 
subsystem. Most notably and due to the existence of the condensation 
phase transition we observe, that volumes of critical droplets flcrit cannot 
be arbitrary small for given subsystem size ^. Their values are bounded 
from below through a ^ dependent bound, which is saturated at the con- 
densation phase transition 

1 



crit - ^ , e'/' (5.29) 



27 mo Acond 

with Acond as given in Eq. ()4.8|) . The volume of critical droplets according 



to Eq. ()5.27p is proportional to and can have any value. At small 
values of h the critical droplet volume is large and the inequality of Eq. 
(j5.29|) is satisfied even for values ^ = L. In this region stochastic tunneling 
takes place with the nucleation of a single droplet. On the other hand, 
at large values of h the critical droplet volume is small and the inequality 
of Eq. ()5.29|) only is satisfied if ^ < L. Consequently tunneling proceeds 
within many subsystems of size ^ i.e., multi-droplet states are encountered. 
For finite size systems with linear extent L and volume we determine 
the magnetic field /id5, where the crossover from single- to multi-droplet 
decay takes place, as 

hDs = 2c2A,ond L-^'^ (5.30) 
and for values h > hos we find the subsystem linear size 
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with C,{h = hos) = L. The interplay of the condensation phase transition 
with the dynamic spinodal was already noted in These authors 
argued however in favor of a logarithmic L - dependency (In L/L)^ instead 
of the dependence of the critical field hps- Recent rigorous argu- 

ments EZIini Qn location of the condensation phase transition exclude 
such leading logarithmic corrections in the finite size behavior of FL{m). 
Our simulation reveals sizable sub-leading finite size corrections for the 
considered lattice sizes, which deserve further studies. With the assump- 
tion, that the multi-droplet physics necessary to derive the location of 
the dynamic spinodal is already contained in the free energy landscape of 
the restricted magnetization equilibrium partition function, logarithmic 
terms should be absent in the leading L dependence of hcs- Further ef- 
forts are necessary to check, whether this assumption holds true or wether 
the dynamical effects invoked in ^ change the L dependence of hns- 

6. CONCLUSION AND OUTLOOK 

The constraint magnetization partition function of the 2D Ising 
model, which we studied at inverse temperature f3 = 0.7, hosts a vari- 
ety of thermodynamic singularities, which qualify as genuine first-order 
phase transitions with free energy barriers, that diverge in the thermody- 
namic limit. The phase space on the torus splits into five disconnected 
sectors: bulk, droplet and strip phases (+/— symmetric states counted 
separately) and, on cube-surface lattices the number of phases is even 
higher, nine in total. On toroidal and cube-surface lattices we observe 
shape transitions, which strongly depend on the choice of the lattice man- 
ifold and, which for finite systems only can be avoided, if the model could 
be formulated on a perfect sphere. Since no such regularization is known, 
we are faced with discontinuous behavior for certain quantities and with 
free energy barriers, which are proportional to the system size L. The 
barrier value and the position of the crystal shape transition on the torus 
(from a droplet to a strip) and the positions and barriers of the "corner 
occupying" droplet shape transitions on the cube-surface, all are very well 
described by a classical droplet description. 

The algorithmic performance of multicanonical ensemble simulations 
suffers from the existence of such shape transitions. Barriers at the discon- 
tinuous phase transitions can not be removed and result into exponential 
slowing down. The earlier general conjecture of random walk behavior for 
multicanonical ensemble simulations in applications to first-order phase 
transitions is falsified for a special case. This interesting algorithmic fact 
untill now has escaped detection, just because free energy barriers in most 
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applications of multicanonical ensemble simulations - though present - 
were not particularly large. There is no doubt, that energy driven first- 
order phase transitions on toroidal lattices with periodic boundaries pos- 
sess similar crystal shape transitions e.g., from a droplet to a strip. In 
higher dimensions than d = 2, e.g. in c? = 3, one faces the additional 
fact, that two phase separation on toroidal boxes in intermediate stages 
also proceeds with the formation of cylinders, which adds additional shape 
transitions to the scenario. We expect, that the performance degradation, 
which is observed in multicanonical ensemble simulations, is a general 
property of broad histogram sampling methods for phase separated sys- 
tems and, that Wang-Landau density of states updating also is affected. 

For toroidal lattices we have obtained a precise finite size parameter- 
ization for droplet and strip free energies - and for the finite size energy 
rounding at the crystal shape transition. Our finding predicts the existence 
of logarithmic phase space factors in the droplet volume, which in addition 
to Gibbs-Thomson and capillary wave fluctuation corrections, contribute 
to the droplet free energy. Such terms reflect the fact, that droplets at 
fixed magnetization may fluctuate to any spatial position and, for the con- 
sidered droplet sizes at the crystal shape transition these corrections are 
actually larger, than e.g. Gibbs-Thomson corrections. A Tolman correc- 
tion to the classical droplet free energy could not be observed, confirming 
the prediction of ^ that the amplitude of this correction vanished for 
systems like the Ising model, which are symmetric under the interchange 
of the two phases. We have conjectured, that translational invariance for 
strip states is broken, because similar logarithmic phase space factors are 
absent for the strip free energy. It would be quite interesting to study the 
dependence of this effect on temperature and dimension. 

The phase separated region of the 2D Ising model is bounded by a 
condensation phase transition, which as we checked at /3 = 0.7 is of dis- 
continuous nature. Within the scope of the present paper we worked out 
the consequences of a simple theoretical model in the one droplet sector, 
which is based on classical droplet theory and Gibbs-Thomson correc- 
tions. Gibbs-Thomson corrected classical droplet theory decomposes the 
free energy into the classical contribution of the droplet, adding fluctua- 
tions through the expansion of the Ginzburg-Landau free energy, which 
in this paper was considered up to second order. This theory performs 
surprisingly well for large droplets. The finite size condensation phase 
transition shift, the gap in the minority droplet size, the ratio of the shift 
over the gap and the nucleation barrier - all seem to approach the second 
order Ginzburg-Landau free energy theoretical prediction. Future stud- 
ies should answer the important theoretical question, whether all of the 
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observed corrections for large, small and smallest droplets can be incor- 
porated through a higher order expansion of the Ginzburg-Landau free 
energy, which as already mentioned requires precise knowledge on the 
shape of the constraint effective potential in vicinity of the bulk. We have 
also presented an argument relating the condensation phase transition to 
the location of the dynamic spinodal, which is of relevance for magnetic 
field switching experiments. 

Thermodynamic singularities associated with the condensation of 
droplets in models with phase separation are expected to exist in any 
dimension - and also for temperature driven first-order phase transitions. 
To our knowledge such transitions have not yet been studied with similar 
methods in numerical simulations and - for most of the time - the effect has 
plainly been overlooked. For purposes of illustration we display in Fig. |^ 
the energy probability distribution function P{E) on a 70^ box from *^ 
in the two-dimensional q = 20 Potts model at the transition temperature. 
The conjectured positions of two asymmetric condensation phase transi- 
tions again is indicated by finite size rounded cusps. The thermodynamic 
properties of these transitions, like their order and their nature in terms 
of droplet and fluctuation degrees have not yet been studied, neither in 
dimension two, nor in higher dimensions. 
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NOTE ADDED IN PROOF 

After completion of this work we became aware of recent work on 
details of the condensation transition by K. Binder (Physica A 319:99 
(2003) , ^cond- mat / 0303651 ) , K. Binder and coworkers (cond- mat/0303642) 
and by M. Biskup, L. Chayes and R. Kotecky ( |cond-mat / 0302373J 
|math.PR70212300 , matirph/0 30203l|. 
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